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Abstract 

We  summarise  research  reported  in  (Efroimsky  2002,  2003;  Efroimsky  &:  Goldreich 
2003a, b)  and  develop  its  application  to  planetary  equations  in  non-inertial  frames. 

Whenever  a  standard  system  of  six  planetary  equations  (in  the  Lagrange,  Delaunay, 
or  other  form)  is  employed,  the  trajectory  resides  on  a  9(N-l)-dimensional  submanifold 
of  the  12(N-l)-dimensional  space  spanned  by  the  orbital  elements  and  their  time  deriva¬ 
tives.  The  freedom  in  choosing  this  submanifold  reveals  an  internal  symmetry  inherent 
in  the  description  of  the  trajectory  by  orbital  elements.  This  freedom  is  analogous  to 
the  gauge  invariance  of  electrodynamics.  In  traditional  derivations  of  the  planetary 
equations  this  freedom  is  removed  by  hand  through  the  introduction  of  the  Lagrange 
constraint,  either  explicitly  (in  the  variation-of-parameters  method)  or  implicitly  (in 
the  Hamilton- Jacobi  approach).  This  constraint  imposes  the  condition  that  the  orbital 
elements  osculate  the  trajectory,  i.e.,  that  both  the  instantaneous  position  and  velocity 
be  fit  by  a  Keplerian  ellipse  (or  hyperbola).  Imposition  of  any  supplementary  con¬ 
straint  different  from  that  of  Lagrange  (but  compatible  with  the  equations  of  motion) 
would  alter  the  mathematical  form  of  the  planetary  equations  without  affecting  the 
physical  trajectory. 

Eor  coordinate-dependent  perturbations,  any  gauge  different  from  that  of  Lagrange 
makes  the  Delaunay  system  non-canonical.  In  a  more  general  case  of  disturbances 
dependent  also  upon  velocities,  there  exists  a  ’’generalised  Lagrange  gauge”  wherein 
the  Delaunay  system  is  symplectic  (and  the  orbital  elements  are  osculating  in  the  phase 
space).  This  gauge  reduces  to  the  regular  Lagrange  gauge  for  perturbations  that  are 
velocity- independent . 

We  provide  a  practical  example  illustrating  how  the  gauge  formalism  considerably 
simplifies  the  calculation  of  satellite  motion  about  an  oblate  processing  planet. 
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Introduction 


1.1  Prefatory  notes 

On  the  6-th  of  November  1766  young  geometer  Joseph-Louis  Lagrange,  ^  invited  from  Turin 
at  d’Alembert’s  recommendation  by  King  Friedrich  the  Second,  succeeded  Euler  as  the  Di¬ 
rector  of  Mathematics  at  the  Berlin  Academy.  Lagrange  held  the  position  for  20  years,  and 
this  fruitful  period  of  his  life  was  marked  by  an  avalanche  of  excellent  results,  and  by  three 
honourable  prizes  received  by  him  from  the  Academie  des  Sciences  of  Paris.  All  three  prizes 
(one  of  which  he  shared  with  Euler)  were  awarded  to  Lagrange  for  his  contributions  to  celes¬ 
tial  mechanics.  Among  these  contributions  was  a  method  initially  developed  by  Lagrange  for 
his  studies  of  planet-perturbed  cometary  orbits  and  only  later  applied  to  planetary  motion 
(Lagrange  1778,  1783,  1788,  1808,  1809,  1810).  The  method  is  based  on  an  elegant  mathe¬ 
matical  novelty,  the  variation  of  parameters  emerging  in  solutions  of  differential  equations. 
The  hrst  known  instances  of  this  tool  being  employed  are  presented  in  a  paper  on  Jupiter’s 
and  Saturn’s  mutual  disturbances,  submitted  by  Euler  to  a  competition  held  by  the  French 
Academy  of  Sciences  (Euler  1748),  and  in  the  treatise  on  the  Lunar  motion,  published  by 
Euler  in  1753  in  St. Petersburg  (Euler  1753).  However  it  was  Lagrange  who  revealed  the  full 
power  of  the  method. 

Below  we  shall  demonstrate  that  the  equations  for  the  instantaneous  orbital  elements 
possess  a  hidden  symmetry  similar  to  the  gauge  symmetry  of  electrodynamics.  Derivation  of 
the  Lagrange  system  involves  a  mathematical  operation  equivalent  to  the  choice  of  a  specihc 
gauge.  As  a  result,  trajectories  get  constrained  to  some  9-dimensional  submanifold  in  the  12- 
dimensional  space  constituted  by  the  orbital  elements  and  their  time  derivatives.  However, 
the  choice  of  this  submanifold  is  essentially  ambiguous,  and  this  ambiguity  gives  birth  to  an 
internal  symmetry.  The  symmetry  is  absent  in  the  2-body  case,  but  comes  into  being  in  the 
N-body  setting  {N  >  3)  where  each  body  follows  an  ellipse  of  varying  shape  whose  time 
evolution  contains  an  inherent  ambiguity. 

For  a  simple  illustration  of  this  point  imagine  two  coplanar  ellipses  sharing  one  focus. 
Suppose  they  rotate  at  different  rates  in  their  common  plane.  Let  a  planet  be  located  at  one  of 
the  intersection  points  of  these  ellipses.  The  values  of  the  elliptic  elements  needed  to  describe 
its  trajectory  would  depend  upon  which  ellipse  was  chosen  to  parameterise  the  orbit.  Either 
set  would  be  equally  legitimate  and  would  faithfully  describe  the  physical  trajectory.  Thus 
we  see  that  there  exists  an  inhnite  number  of  ways  of  dividing  the  actual  planet’s  movement 
into  motion  along  its  orbit  and  the  simultaneous  evolution  of  the  orbit.  Although  the  physical 
trajectory  is  unique,  its  description  (parametrisation  in  terms  of  Kepler’s  elements)  is  not. 
A  map  between  two  different  (though  physically-equivalent)  sets  of  orbital  elements  is  a 
symmetry  transformation  (a  gauge  transformation,  in  physicists’  jargon). 

Lagrange  never  dwelled  on  that  point.  However,  in  his  treatment  (based  on  direct  ap¬ 
plication  of  the  method  of  variation  of  parameters)  he  passingly  introduced  a  convenient 
mathematical  condition  which  removed  the  said  ambiguity.  This  condition  and  possible 
alternatives  to  it  will  be  the  topics  of  Sections  1  -  3  of  this  paper. 

^The  real  name  of  the  young  man  invited  in  1766  to  the  Prussian  court  was  Giuseppe  Lodovico  Lagrangia. 
It  was  only  several  years  down  the  road  that  he  became  Joseph-Louis  Lagrange. 
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In  1834  -  1835  Hamilton  put  forward  his  theory  of  canonical  transformations.  Several 
years  later  this  approach  was  furthered  by  Jacobi  who  brought  Hamilton’s  technique  into 
astronomy  and,  thereby,  worked  out  a  new  method  of  deriving  the  planetary  equations 
(Jacobi  1866),  a  method  that  was  soon  accepted  as  standard.  Though  the  mathematical 
content  of  the  Hamilton-Jacobi  theory  is  impeccably  correct,  its  application  to  astronomy 
contains  a  long  overlooked  aspect  that  needs  attention.  This  aspect  is:  where  is  the  Lagrange 
constraint  tacitly  imposed,  and  what  happens  if  we  impose  a  different  constraint?  This  issue 
will  be  addressed  in  Section  4  of  our  paper. 

1.2  Osculating  Elements  vs  Orbital  Elements 

We  start  in  the  spirit  of  Lagrange.  Before  addressing  the  N-particle  case,  Lagrange  referred 
to  the  reduced  2-body  problem. 


r  -h 


/i  r 


=  0 


(1) 


r  = 


'^planet 


h  —  ^(j^planet  T  ^^sun) 


whose  generic  solution,  a  Keplerian  ellipse  or  a  hyperbola,  can  be  expressed,  in  some  fixed 
Cartesian  frame,  as 


where 


=  g(Cl,--,C'6,  t) 


i 


C=const 


(2) 

(3) 


Since  the  problem  ((T))  is  constituted  by  three  second-order  differential  equations,  its  general 
solution  naturally  contains  six  adjustable  constants  Cj.  At  this  point  it  is  irrelevant  which 
particular  set  of  the  adjustable  parameters  is  employed.  (It  may  be,  for  example,  a  set 
of  Lagrange  or  Delaunay  orbital  elements  or,  alternatively,  a  set  of  initial  values  of  the 
coordinates  and  velocities.) 

Following  Lagrange  (1808,  1809,  1810),  we  employ  /  as  an  ansatz  for  a  solution  of  the 
N-particle  problem,  the  disturbing  force  acting  at  a  particle  being  denoted  by  AF 


r  -h 


/i  r 


AF 


Now  the  ’’constants”  become  time  dependent: 

r  =  m(t),...,C6(t),t) 


(4) 

(5) 


whence  the  velocity 

dt  dt  ^  dCi  dt  ^  i  9Ci  dt  ’ 

^Our  treatment  covers  disturbing  forces  AF(f,  r,  t)  that  are  arbitrary  vector- valued  functions  of  position, 
velocity,  and  time. 
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acquires  a  new  input,  Y.{df/dCi){dCi/dt)  . 

Substitution  of  ©  and  into  the  perturbed  equation  of  motion  0  leads  to  three  inde¬ 
pendent  scalar  second-order  differential  equations  which  contain  one  independent  parameter, 
time,  and  six  functions  Ci{t)  .  These  are  to  be  found  from  the  said  three  equations,  and 
this  makes  the  essence  of  the  variation-of-parameters  (VOP)  method.  However,  the  latter 
task  cannot  be  accomplished  in  a  unique  way  because  the  number  of  variables  exceeds,  by 
three,  the  number  of  equations.  Thence,  though  the  physical  trajectory  (comprised  by  the 
locus  of  points  in  the  Cartesian  frame  and  by  the  values  of  velocity  at  each  of  these  points) 
is  unique,  its  parametrisation  through  the  orbital  elements  is  ambiguous.  This  circumstance 
was  appreciated  by  Lagrange  who  amended  the  system,  by  hand,  with  three  independent 
conditions, 

dCi  dt 

and  thus  made  it  solvable.  His  choice  of  constraints  was  motivated  by  both  physical  consider¬ 
ations  and  mathematical  expediency.  Since,  physically,  the  set  of  functions  (C'i(f),  ...,0^(1)) 
can  be  interpreted  as  parameters  of  an  instantaneous  ellipse,  in  a  bound-orbit  case,  or  an 
instantaneous  hyperbola,  in  a  fly-by  situation,  Lagrange  found  it  natural  to  make  the  instan¬ 
taneous  orbital  elements  Ci  osculating.  His  constraint  m  fixes  the  instantaneous  ellipse 
(or  hyperbola)  defined  by  (Ci(f),  ...,0^(1))  such  that,  at  each  moment  of  time,  it  coincides 
with  the  unperturbed  (two-body)  orbit  that  the  body  would  follow  if  the  disturbances  were 
to  cease  instantaneously.  This  way,  Lagrange  restricted  his  use  of  the  orbital  elements  to 
elements  that  osculate  in  the  reference  frame  wherein  ansatz  0  is  employed.  Lagrange  never 
explored  alternative  options;  he  simply  imposed  o  and  used  it  to  derive  his  famous  system 
of  equations  for  orbital  elements. 

Such  a  restriction,  though  physically  motivated,  is  completely  arbitrary  from  the  math¬ 
ematical  point  of  view.  While  the  imposition  of  m  considerably  simplihes  the  subsequent 
calculations  it  in  no  way  influences  the  shape  of  the  physical  trajectory  and  the  rate  of 
motion  along  it.  A  choice  of  any  other  supplementary  constraint 


E 


df  da 


dCi  dt 


=  ^(C'wW) 


(8) 


$  being  an  arbitrary  function  of  time  and  parameters  Ci,  would  lead  to  the  same  physical 
orbit  and  the  same  velocities.^  Substitution  of  the  Lagrange  constraint  m  by  its  general¬ 
isation  (0  would  not  influence  the  motion  of  the  body  but  would  alter  its  mathematical 
description  (i.e.,  would  entail  different  solutions  for  the  orbital  elements).  Such  invariance 
of  a  physical  picture  under  a  change  of  parametrisation  is  called  gauge  freedom  or  gauge 
symmetry.  It  parallels  a  similar  phenomenon  well  known  in  electrodynamics  and,  therefore, 
has  similar  consequences.  On  the  one  hand,  a  good  choice  of  gauge  often  simplihes  solution 
of  the  equations  of  motion.  (In  Section  3  we  provide  a  specihc  application  to  motion  in  an 
accelerated  coordinate  system.)  On  the  other  hand,  a  so-called  “gauge  shift”  will  occur  in 
the  course  of  numerical  orbit  computation  (Murison  &  Efroimsky  2003). 

Derivation  of  the  conventional  Lagrange  and  Delaunay  planetary  equations  by  the  VOP 
method  incorporates  the  Lagrange  constraint  (Brouwer  &  Clemence  1961).  Both  systems 

^In  principle,  one  may  endow  $  also  with  dependence  upon  the  parameters’  time  derivatives  of  all  orders. 
This  would  yield  higher-than-first-order  time  derivatives  of  the  Ci  in  subsequent  developments  requiring 
additional  initial  conditions,  beyond  those  on  r  and  r,  to  be  specified  in  order  to  close  the  system.  We  avoid 
this  unnecessary  complication  by  restricting  $  to  be  a  function  of  time  and  the  Ci  . 
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of  equations  get  altered  under  a  different  gauge  choice  as  we  now  show.  The  essence  of  a 
derivation  suitable  for  a  general  choice  of  gauge  starts  with  (jH))  from  which  the  formula  for 
the  acceleration  follows: 


d'^r 

dt‘^ 


^  _L  V  —  ^  I  ^ 

dt  ^  dCi  dt  dt 

I 


dH  ^  dCi  dt  dt 


Together  with  the  equation  of  motion  it  yields: 

_L  ii  /  +  V  —  ^  ^  ^ 

dt'^  r  ^  dCi  dt  dt 

2  ^ 


(9) 


(10) 


The  vector  function  /  was  from  the  beginning  introduced  as  a  Keplerian  solution  to  the 
two-body  problem;  hence  it  must  obey  the  unperturbed  equation  O-  On  these  grounds  the 
above  formula  simplihes  to: 


^  dj  dCj  _  d^ 

^  dCi  dt  dt 


(11) 


This  equation  describes  the  perturbed  motion  in  terms  of  the  orbital  elements.  Together 
with  constraint  (jHl)  it  constitutes  a  well-dehned  system  which  may  be  solved  with  respect 
to  dCi/dt  .  An  easy  way  to  do  this  is  to  use  the  elegant  trick  suggested  by  Lagrange:  to 
multiply  the  equation  of  motion  by  df  /dCn  and  to  multiply  the  constraint  by  —  dg/dCn  ■ 
The  former  operation  results  in 


^  ^  ^  AF  _  ^  ^ 

dCn  \~  dCj  dt  j  dCn  dCn  dt 

while  the  latter  gives 

_  dg  ^ 

dCn  dCj  dt  )  dCn 

Having  summed  these  two  equalities,  we  arrive  to: 


i:  ic„  Qi 


dCj 

dt 


dCn  dCn  dt  dCn 


(12) 


(13) 


(14) 


|C„  c,] 

brackets: 


standing  for  the  unperturbed  (i.e.,  dehned  as  in  the  two-body  case)  Lagrange 


^  ~  dCn  dCj  dCj  dCn 


(15) 


It  was  agreed  above  that  $  is  a  function  of  time  and  of  the  ’’constants”  Ci  ,  but  not 
of  their  time  derivatives.  Under  this  convention,  the  above  equation  may  be  shaped  into  a 
more  convenient  form: 


E 


[a  c,] 


df  d$ 


dt 


df  ^  ^  g 

dCn  dCn  dt  dCn 


(16) 
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This  is  the  most  general  form  of  the  gauge-invariant  perturbation  equations  of  celestial 
mechanics.  In  the  Lagrange  gauge,  when  the  $  terms  are  absent,  we  can  obtain  an 
immediate  solution  for  the  individual  dCi/dt  by  exploiting  the  well  known  expression  for 
the  Poisson-bracket  matrix  which  is  inverse  to  the  Lagrange-bracket  one  and  is  offered  in 
the  literature  for  the  two-body  problem.  (Be  mindful  that  our  brackets  (jl5|l  are  defined  in 
the  same  manner  as  in  the  two-body  case;  they  contain  only  functions  f  and  g  that  are 
defined  in  the  unperturbed,  two-body,  setting.)  In  an  arbitrary  gauge  the  presence  of  the 
term  proportional  to  d^/dCj  on  the  left-hand  side  of  (ll(i|l  complicates  the  solution  for 
dCi/dt,  but  only  to  the  extent  of  requiring  the  resolution  of  a  set  of  six  simultaneous  linear 
algebraic  equations. 

To  draw  to  a  close,  we  again  emphasise  that,  for  fixed  interactions  and  initial  conditions, 
all  possible  (i.e.,  compatible  and  sufficient)  choices  of  gauge  conditions  expressed  by  the 
vector  function  $  lead  to  a  physically  equivalent  picture.  In  other  words,  the  real  trajectory 
is  invariant  under  reparametrisations  permitted  by  the  ambiguity  of  the  choice  of  gauge.  This 
invariance  has  the  following  meaning.  Suppose  the  equations  of  motion  for  ,  with 

some  gauge  condition  $  imposed,  render  the  solution  .  The  same  equations  of 

motion,  with  another  gauge  $  enforced,  furnish  the  solution  Ci{t)  that  has  a  different 
functional  form.  Despite  this  difference,  both  solutions,  Ci{t)  and  Ci{t)  ,  when  substituted 
back  in  yield  the  same  curve  r(t)  with  the  same  velocities  r(f)  .  In  mathematics  this 
situation  is  called  a  fibre  bundle,  and  it  gives  birth  to  a  1-to-l  map  of  Ci{t)  onto  Ci{t)  , 
which  is  merely  a  reparametrisation.  In  physics  this  map  is  called  a  gauge  transformation. 
The  entire  set  of  these  reparametrisations  constitute  a  group  of  symmetry  known  as  a  gauge 
group. 

Just  as  in  electrodynamics,  where  the  fields  E  and  B  stay  invariant  under  gradient  trans¬ 
formations  of  the  4-potential  ,  so  the  invariance  of  the  orbit  implements  itself  through 
the  form-invariance  of  expression  ©  under  the  afore  mentioned  map.  The  vector  r  and  its 
full  time  derivative  r  ,  play  the  role  of  the  physical  fields  E  and  B ,  while  the  Keplerian 
coordinates  Ci, ...,  Cg  play  the  role  of  the  4-potential  . 

A  comprehensive  discussion  of  all  the  above-raised  issues  can  be  found  in  Efroimsky 
(2002,  2003).  Interconnection  between  the  internal  symmetry  and  multiple  time  scales  in 
celestial  mechanics  is  touched  upon  in  Newman  &  Efroimsky  (2003). 


2  Planetary  equations  in  an  arbitrary  gauge 


2.1  Lagrangian  and  Hamiltonian  Perturbations 


We  can  proceed  further  by  restricting  the  class  of  perturbations  we  consider  to  those  in  which 
AF  is  derivable  from  a  perturbed  Lagrangian.  This  restricted  class  is  still  broad  enough  to 
encompass  most  applications  of  celestial-mechanics  perturbation  theory. 

Let  the  unperturbed  dynamics  be  described  by  Lagrangian  Lo(r,  r,  t)  =  r  /2  — 
f/(r  ,  t),  canonical  momentum  p  =  r  and  Hamiltonian  Ho{r,  p,  f)  =  p  ^/2  -|-  f/(r  ,  t)  . 
Disturbed  motion  will  be  described  by  the  new,  perturbed,  functions: 

2 

L(r,  r,  t)  =  Lo  +  AL  =  ^ - [/(r,  t)  +  AL(r,  r,  t)  ,  (17) 


P 


dAL 

r  +  - ^ 

dr 


(18) 
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p 


H{r,  p,  t)  =  pr  -  L  =  — — \-  U  +  AH  ,  Aif(r,  p,  t)  =  -  AL  - 


fdALV 
V  (9r  / 


(19) 


AH  being  introduced  as  a  variation  of  the  functional  form,  i.e.,  as  AH  =  H{r,  p,  t)  — 
Ho{r,  p,  t)  .  The  Euler-Lagrange  equations  written  for  the  perturbed  Lagrangian  (fTTjl  will 
give: 


r  =  — 


dU 

dr 


+  AF 


(20) 


where  the  new  term 


dAL  _  d  fdAL\ 
a?  dt  {  gf  J 


(21) 


is  the  disturbing  force.  Expression  m  reveals  that  whenever  the  Lagrangian  perturbation 
is  velocity  independent,  it  plays  the  role  of  the  disturbing  function.  Generally,  though,  the 
disturbing  force  is  not  equal  to  the  gradient  of  AL  but  has  an  extra  term  generated  by  the 
velocity  dependence. 

Examples  in  which  a  velocity  dependent  AL  has  been  used  in  a  celestial  mechanics 
setting  include:  the  treatment  of  inertial  forces  in  a  coordinate  system  tied  to  the  spin  axis 
of  a  processing  planet  (Goldreich  1965);  the  velocity  dependent  corrections  to  Newton’s  law 
of  gravity  in  the  relativistic  two-body  problem  (Brumberg  1992). 


2.2  Gauge- invariant  planetary  equations 


When  the  expression  m  for  the  most  general  force  emerging  within  the  Lagrangian  for¬ 
malism  is  substituted  into  the  gauge-invariant  perturbation  equation  (d,  it  yields: 


y  ‘  dt  sc„  dr  dc„  dt[  ^  gf  j  gc„ 

Since,  for  a  velocity-dependent  disturbance,  the  chain  rule 

dAL  dAL  df  dAL  dr  dAL  df  dAL  d{g  +  ^>) 

takes  place,  formula  (1^  may  be  re-shaped  into 

y  "  dt  gCn  gp  dc„  dt  gc„j  t  9f  / 


(22) 


(23) 


(24) 


• 

Next  we  group  terms  so  that  the  gauge  function  $  everywhere  appears  added  to  d{AL)  / dr  , 
and  bring  the  term  proportional  to  dCj / dt  to  the  left  hand  side  of  the  equation: 


E 


( 


[a 


df 

dCn  dCj 


/dAL 
\  dr 


dt 


d 


AL  + 


m' 


dg  df  d  dAL  ^  \  d AL\ 

dCn  dCn  dt  dr  dCn  y  \  dr  j 


(25) 
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These  modifications  help  us  to  recognise  the  special  nature  of  the  gauge  $  =  —  d{AL) / dr 
which  will  be  the  subject  of  discussion  in  the  next  subsection. 

Contrast  (123)  with  (USD:  while  (fT3  expresses  the  VOP  method  in  the  most  general  form 
it  can  have  in  terms  of  the  disturbing  force  AF(r,  r,  t)  ,  equation  (I25|l  renders  the  most 
general  form  in  the  language  of  a  Lagrangian  perturbation  AL(r,  r,  t). 

The  applicability  of  so  generalised  planetary  equations  in  analytical  calculations  is  com¬ 
plicated  by  the  nontrivial  nature  of  the  left-hand  sides  of  (16)  and  (25).  Nevertheless,  the 
structure  of  these  left-hand  sides  leaves  room  for  analytical  simplification  in  particular  situ¬ 
ations.  One  such  situation  is  when  the  gauge  is  chosen  to  be 

-  d  AL  ,  ,  , 

$  = - ^  +  r}{r,  t)  ,  (26) 

dr 

ff{r,  t)  being  an  arbitrary  vector  function  linear  in  r.  (It  may  be,  for  example,  proportional 
to  r  or  may  be  equal,  say,  to  a  cross  product  of  r  by  some  time-dependent  vector.)  Under 
these  circumstances  the  left-hand  side  in  (25)  reduces  to  the  Lagrange  brackets.  The  situation 
becomes  especially  simple  when  dAL / dr  happens  to  be  linear  in  r  ,  in  which  case  we  may 
put  ff{r,  t)  =  dAL /dr  and,  thus,  employ  the  trivial  Lagrange  gauge  $  =  0  instead  of 
the  generalised  Lagrange  gauge.  We  shall  encounter  one  such  example  in  section  3.4. 

As  already  stressed  above,  the  Lagrange  brackets  are  gauge-invariant  because  functions 
f  and  g  used  in  are  defined  within  the  unperturbed,  two-body,  problem  (P-El)  that  lacks 
gauge  freedom.  For  this  reason,  one  may  exploit,  to  solve  (123,  the  well-known  expression 
for  the  inverse  of  matrix  [Ci  Cj\  .  Its  elements  are  simplest  (and  are  either  zero  or  unity) 
when  one  chooses  as  the  ’’constants”  the  Delaunay  set  of  orbital  variables  (Plummer  1918): 


a  =  {  L  ,  G ,  if ,  M,  ,  0^  ,  D  } 


L  =  ^/Jla  , 


G  =  \jfia  (1  —  e^) 


e^)  cos  i 


(27) 


where  /x  =  G  {m  sun  +  ^planet)  and  (e,  a,  Mq  ,  a; ,  D  ,  i)  are  the  Keplerian  elements:  e 
and  a  are  the  eccentricity  and  major  semiaxis.  Mo  is  the  mean  anomaly  at  epoch,  and  the 
Euler  angles  u  ,  Q ,  i  are  the  the  argument  of  pericentre,  the  longitude  of  the  ascending 
node,  and  the  inclination,  respectively. 

The  simple  forms  of  the  Lagrange  and  Poisson  brackets  in  Delaunay  elements  is  the  proof 
of  these  elements’  canonicity  in  the  unperturbed,  two-body,  problem:  the  Delaunay  elements 
give  birth  to  three  canonical  pairs  {Qi,  Pi)  corresponding  to  a  vanishing  Hamiltonian: 
(L  ,  —  Mo) ,  (G  ,  —  a;) ,  {H  ,  —Q).  In  a  perturbed  setting,  when  only  a  position-dependent 
disturbing  function  R{r,  t)  is  ’’turned  on”,  it  can  be  expressed  through  the  Lagrangian  and 
Hamiltonian  perturbations  in  a  simple  manner,  R{r,  t)  =  AL{r,  t)  =  —AH{r,  t)  ,  as  can 
be  seen  from  the  formulae  presented  in  the  previous  subsection.  Under  these  circumstances, 
the  Delaunay  elements  remain  canonical,  provided  the  Lagrange  gauge  is  imposed  (Brouwer 
&  Clemence  1961).  This  long  known  fact  can  also  be  derived  from  our  equation  (I25|l:  if  we 
put  $  =  0  and  assume  AL  velocity-independent,  we  arrive  to 


where 


i:  ic„  Q] 


dt 


dAL 


(28) 


AL  =  AL{f(C,t),t)  = 


AH  {f(C,t),  t) 


(29) 


This,  in  its  turn,  results  in  the  well  known  Lagrange  system  of  planetary  equations,  provided 
the  parameters  Ci  are  chosen  as  the  Kepler  elements.  In  case  they  are  chosen  as  the 
Delaunay  elements,  then  leads  to  the  standard  Delaunay  equations,  i.e.,  to  a  symplectic 
system  wherein  the  pairs  {L  ,  —  Mo) ,  (G  ,  —  a;) ,  {H  ,  —  D)  again  play  the  role  of  canonical 
variables,  but  the  Hamiltonian,  in  distinction  to  the  unperturbed  case,  no  longer  vanishes, 
instead  being  equal  to  AiL  =  —  AL  . 

In  our  more  general  case,  the  perturbation  depends  also  upon  velocities  (and,  therefore, 
AL  is  no  longer  equal  to  —AH  ).  Beside  this,  the  gauge  $  is  set  arbitrary.  As 
demonstrated  in  Efroimsky  (2002,  2003),  under  these  circumstances  the  gauge- invariant 
Delaunay-type  system  is  no  longer  symplectic.  However,  it  turns  out  that  this  system  regains 
the  canonical  form  in  one  special  gauge,  one  that  coincides  with  the  Lagrange  gauge  when  the 
perturbation  bears  no  velocity  dependence.  The  issue  is  explained  at  length  in  our  previous 
papers  (Efroimsky  &  Goldreich  2003a,b).  Here  we  offer  a  brief  synopsis  of  this  study. 


2.3  Generalised  Lagrange  gauge  wherein  the  Delaunay-type 
system  is  canonical 


Equation  (0^  was  cast  in  the  shape  of  (EiD  not  only  to  demonstrate  the  special  nature  of 
the  gauge 


$ 


dAL 

dr 


(30) 


but  also  to  single  out  the  terms  in  the  square  brackets  on  the  right-hand  side  of  together, 
these  terms  give  exactly  the  Hamiltonian  perturbation.  Thus  we  come  to  the  conclusion  that 
in  the  special  gauge  dSOD  our  equation  (ESI)  simplifies  to 


i:  [Cn  C,] 


dCj 

dt 


d  AH 
dCn 


(31) 


As  emphasised  in  the  preceding  subsection,  the  gauge  invariance  of  definition  ESD  enables  us 
to  use  the  standard  (Lagrange-gauge)  expressions  for  [G„  Cj\~^  to  get  the  planetary  equations 
from  (jSH). 

Comparing  m  with  (EHl),  we  see  that  in  the  general  case  of  an  arbitrary  AL(r,  r,  t)  one 
arrives  from  EH)  to  the  same  equations  as  from  ()28|1.  except  that  now  they  contain  —  AH 
instead  of  AL.  When  the  orbit  is  parametrised  by  the  Delaunay  variables,  those  equations 
take  the  form: 


dL 

dAH 

d{  -  Mo) 

dAH 

(32) 

dt 

d{-Mo)  ’ 

dt 

dL 

dG 

dAH 

d{  —  oj) 

dAH 

(33) 

dt 

d{-u;) 

’  dt 

dG 

dH 

dAH 

d{-Q) 

dAH 

(34) 

dt 

d{-n)  ’ 

dt 

dH 

which  is  a  symplectic  system.^  For  this  reason  we  name  this  special  gauge  the  ’’generalised 
Lagrange  gauge”.  In  any  different  gauge  $  the  equations  for  the  Delaunay  variables 

^  In  this  system  H  stands  not  for  the  Hamiltonian  but  for  one  of  the  Delaunay  elements. 
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would  contain  ^-dependent  terms  and  would  not  be  symplectic.  (Those  gauge-invariant 
equations,  for  both  Lagrange  and  Delaunay  elements  are  presented  in  Efroimsky  (2002,  2003) 
and  Efroimsky  &  Goldreich  (2003a,b).).  This  analysis  proves  the  following  THEOREM : 
Though  the  gauge-invariant  eqnations  for  Delaunay  elements  are,  generally,  not 
canonical,  they  become  canonical  in  the  ’’generalised  Lagrange  gauge”.  That  this 
Theorem  is  not  merely  a  mathematical  coincidence  but  has  deep  reasons  beneath  it  will  be 
shown  in  Section  4  where  the  subject  is  approached  from  the  Hamilton- Jacobi  viewpoint. 

The  above  Theorem  gives  one  example  of  the  gauge  formalism  being  of  use:  an  appro¬ 
priate  choice  of  gauge  can  considerably  simplify  the  planetary  equations  (in  this  particular 
case,  it  makes  them  canonical). 

According  to  (UHl),  the  momentum  can  be  written  as 


P 


^  dAL 

r 

dr 


g  -h  $  -h 


dAL 

dr 


which,  in  the  generalised  Lagrange  gauge  O,  simply  reduces  to 


P  =  g 


(35) 


(36) 


Vector  g  was  introduced  back  in  (0  -  El)  to  denote  the  functional  dependence  of  the 
unperturbed  velocity  upon  the  time  and  the  parameters  Ci  .  In  the  unperturbed,  two-body, 
setting  this  velocity  is  equal  to  the  momentum  canonically  conjugate  to  the  position  r  (this 
is  obvious  from  eg,  for  zero  AL  ).  This  way,  in  the  unperturbed  case  equality  (Ei  is 
fulhlled  trivially.  The  fact  that  it  remains  valid  also  under  perturbation  means  that,  in  the 
said  gauge,  the  canonical  momentum  in  the  disturbed  setting  is  the  same  function  of  time 
and  ’’constants”  as  in  the  unperturbed,  two-body,  case.  Thus  we  have  established  that  the 
instantaneous  Keplerian  ellipses  (hyperbolae)  defined  in  gauge  En|l  osculate  the  trajectory 
in  phase  space. 

Not  surprisingly,  the  generalised  Lagrange  gauge  (EDI)  reduces  to  =  0  in  the  simple 
case  of  velocity-independent  disturbances. 


3  Gauge  Freedom  and  Freedom  of  Frame  Choice 

3.1  Osculating  ellipses  described  in  different  frames  of  reference. 

The  essence  of  the  VOP  method  in  celestial  mechanics  is  the  following.  A  generic  two-body- 
problem  solution  expressed  by 


r 

= 

fic,t)  , 

(37) 

^df\ 

V*) 

= 

g  {C,  t)  , 

(38) 

c 

/  dg'^ 
[dtj 

Iv 

h  / 

p  f  ’ 

(39) 

is  employed  as  an  ansatz  to  solve  the  disturbed  problem: 

f  =  f{C{t),  t)  , 


(40) 
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r 


dt  dCi  dt 


g  +  ^ 


(41) 


"  di  ,  dg  dCi  d^ 
r  =  +  ttt:; - 7-  +  -rr 


dt  dCi  dt  dt 


(42) 


P  f  dCi  dt  dt 


fif  dgdCi  d^ 

n  '  r\y^  1,  '  1, 


As  evident  from  (EH),  our  choice  of  a  particular  gauge  is  equivalent  to  decomposition  of 
the  physical  motion  into  a  movement  with  velocity  g  along  the  instantaneous  ellipse  (or 
hyperbola,  in  the  fly-by  case),  and  a  movement  associated  with  the  ellipse’s  (or  hyperbola’s) 
deformation  that  goes  at  the  rate  .  It  is  then  tempting  to  state  that  a  choice  of  gauge  is 
equivalent  to  a  choice  of  an  instantaneous  comoving  reference  frame  wherein  to  describe  the 
motion.  Such  an  interpretation  is,  however,  incomplete.  Beside  the  fact  that  we  decouple 
the  physical  velocity  in  a  certain  proportion  between  g  and  ,  it  also  matters  which 
physical  velocity  (i.e.,  velocity  relative  to  what  frame)  is  decoupled  in  this  proportion.  In 
other  words,  our  choice  of  the  gauge  does  not  yet  exhaust  all  freedom:  we  can  still  choose  in 
what  frame  to  write  ansatz  PH).  We  may  write  it  in  inertial  axes  or  in  some  accelerated 
system.  For  example,  in  the  case  of  a  satellite  orbiting  an  accelerated  and  precessing  planet 
it  is  convenient  to  write  the  ansatz  for  the  planet-related  position  vector. 

The  above  kinematic  formulae  PH)  -  dH  do  not  yet  contain  information  about  our  choice 
of  the  reference  system  in  which  we  implement  the  VOP  method.  This  information  shows 
up  at  the  next  stage,  when  expression  EH  is  inserted  into  the  dynamical  equation  of  motion 
r  =  —  (pr/r^)  -|-  AF  to  yield 


dg  dCi  ^  d^ 


dAL  _  d  fdAL 
dr  dt  I 


(43) 


dCi  dt  dt 


Complete  information  about  the  reference  system  in  which  we  put  the  VOP  method  to  work 
(and,  therefore,  in  which  we  dehne  the  orbital  elements  Ci )  is  contained  in  the  expression 
for  the  perturbation  force  AF  .  For  example,  if  the  operation  is  carried  out  in  an  inertial 
coordinate  system,  AF  contains  physical  forces  solely.  However,  if  we  wish  to  implement 
the  VOP  approach  in  a  frame  moving  with  a  linear  acceleration  a  ,  then  AF  also  contains 
the  inertial  force  —  a  .  In  case  this  coordinate  system  rotates  relative  to  inertial  ones  at  a 
rate  fl ,  then  AF  also  includes  the  inertial  terms  — 2/2x  r  —  /Xx  r  —  fl  x  {fl  x  r)  .  In 
considering  the  motion  of  a  satellite  orbiting  an  oblate  precessing  planet  it  is  most  reasonable, 
though  not  obligatory,  to  apply  the  method  (i.e.,  to  dehne  the  time  derivative)  in  axes  that 
precess  with  the  planet.  However,  this  reasonable  choice  of  coordinate  system  still  leaves  us 
with  the  freedom  of  gauge  nomination. 

3.2  Relevant  Example 

Gauge  freedom  of  the  perturbation  equations  of  celestial  mechanics  hnds  an  immediate 
practical  implementation  in  the  description  of  test  particle  motion  around  an  precessing 
oblate  planet  (Goldreich,  1965).  It  is  trivial  to  extend  this  to  account  for  acceleration  of  the 
planet’s  centre  of  mass. 
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Our  starting  point  is  the  equation  of  motion  in  the  inertial  frame 


r  = 


dV 

(9r 


(44) 


where  U  is  the  total  gravitational  potential  and  time  derivatives  in  the  inertial  axes  are 
denoted  by  primes.  Suppose  that  the  planet’s  spin  axis  precesses  at  angular  rate  /A(f)  and 
that  the  acceleration  of  its  centre  of  mass  is  given  by  a(f) .  In  a  coordinate  system  attached 
to  the  planet’s  centre  of  mass  and  processing  with  it,  inertial  forces  modify  the  equation  of 
motion  so  that  it  assumes  the  form: 

"  dU  ^  u  ^  ^  ^  ^ 

r  =  —  2/2  xr  —  /2xr  —  /2x(/2xr)  —  a,  (45) 

ar 

time  derivatives  in  the  accelerated  frame  being  denoted  by  dots. 

To  implement  the  VOP  approach  in  terms  of  the  orbital  elements  defined  in  the  accel¬ 
erated  frame,  we  note  that  the  disturbing  force  on  the  right-hand  side  of  (gSl)  is  generated 
according  to  (ED  by 


AL  (^r,  r,  =  R  +  r-(/2  x  r)  -|-  ^  (/2  x  r)-(/2  x  r)  —  a  ■  r  , 


(46) 


where  we  denote  by  R{r  ,  t)  the  gravitational-potential  perturbation  (which  the  perturbation 
of  the  overall  gravitational  potential  U).  Since 


dAL 

dr 


=  H  X  r 


the  corresponding  Hamiltonian  perturbation  reads: 

AH  =  - 


ifdAlV 

AL  +  -  ‘  ' 


2  \  dr  J 


(47) 


(48) 


' R  +  p-  {fix  r) 


a  ■  r  = 


i?-F(rxp)-/2-a-r]  , 


with  vector  J  =  r  x  p  being  the  satellite’s  orbital  angular  momentum  in  the  inertial 
frame. 

According  to  (1^  and  (HZD,  the  momentum  can  be  written  as 


p  =  g  +  $  +  /2x/. 


whence  the  Hamiltonian  perturbation  becomes 


AH 


(49) 


(50) 


3.3  Elements  defined  in  an  accelerated,  rotating  frame 
that  osculate  in  the  comoving  inertial  frame 

In  this  subsection  we  recall  a  calculation  carried  out  by  Goldreich  (1965)  and  Brumberg, 
Evdokimova  &  Kochina  (1971)  and  demonstrate  that  it  may  be  interpreted  as  an  example 
of  nontrivial  gauge  hxing. 
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Let  us  implement  the  VOP  method  in  a  frame  that  is  accelerating  at  rate  a  and  rotating 
at  angular  rate  p.  relative  to  some  inertial  system  S.  This  means  that,  in  the  VOP  equation 
(m,  AL  is  given  by  formula  (HH|)  and  AH  by  dsni). 

We  now  choose  to  describe  the  motion  in  the  generalised  Lagrange  gauge  dSHD,  so  the 
expression  +  //  x  r ^  on  the  right-hand  side  of  (j5()|l  vanishes  (as  follows  from  (HZD), 
and  the  expression  for  AH  in  terms  of  /  and  g  has  the  form: 


AH 


R(f,  «)  +  M'  (/ X  g)  -  5  -  / 


(51) 


At  the  same  time,  the  generic  expression  for  the  VOP  given  by  (I25|l  simplihes  to  (ED). 
Insertion  of  m  therein  leads  us  to 


[Cr  a] 


dt 


d 

Wr 


R(f,  i)  +  (/ X  g)  -  a  -  / 


(52) 


Interestingly,  this  equation  does  not  contain  /X  even  though  it  is  valid  for  non-uniform 
precession. 

As  explained  in  subsection  2.3,  in  the  generalised  Lagrange  gauge  the  vector  g  is  equal 
to  the  canonical  momentum  p  =  r  -|-  dAL/dr.  In  the  case  when  the  velocity  dependence 
of  AL  is  called  into  being  by  inertial  forces,  the  momentum  is,  according  to  dUD, 


^  U  LSLj  ^ 

p  =  r  -|-  - ^  =  r  -|-  p  X  r  ,  (53) 

dr 

which  is  the  particle’s  velocity  relative  to  the  inertial  frame  comoving  with  the  accelerated, 
rotating  frame.  In  this  sense  we  may  say  that  our  elements  are  dehned  in  the  accelerated, 
rotating  frame  but  osculate  in  the  comoving  inertial  one. 

In  the  appendix  we  provide  explicit  expression  for  each  of  the  partial  derivatives  of  p  -  J 
that  appears  in  the  planetary  equations  (IS21). 


3.4  Elements  defined  in  the  accelerated,  rotating  frame, 
that  osculate  in  this  frame 


Here  we  not  only  dehne  the  elements  in  the  accelerated,  rotating  frame,  but  we  also  make 
them  osculate  in  this  system,  i.e.,  we  make  them  satisfy  =  0.  In  this  gauge,  expression 
(EOD  takes  the  following  form: 


AH  =  - 


R{f,  f)  -f  p  ■  (/ X  g)  +  (p  X  /)  ■  (p  X  /)  -  a  ■  /] 

while  equation  (I2ni),  after  some  algebra,®  looks  like  this: 

dQ  dAH 


dC„ 


X  i  -  f  X 


dC„ 


[Cn  a 


P-  \  f  X 


dt 


dC„ 


dCr, 


d 


(54) 


(65) 


®  Due  to  63),  the  second  term  on  the  left-hand  side  in  (I25II  is  proportional  to  [d(/2  x  f)/dCj]  Cj  =  px  $ 
and,  therefore,  vanishes.  The  second  term  on  the  right-hand  side  simplifies  in  accordance  with  the  simple 
rule  A  •  (B  X  C)  =  (C  x  A)  •  B  .  (See  Efroimsky  &  (Goldreich  2003b).) 
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When  substituting  dMl)  into  it  is  convenient  to  rent  the  expression  for  AH  apart  and 
to  group  the  term  {p,  x  /)  ■  (m  X  /)  with  the  last  term  on  the  right-hand  side  of 


[Cn  a] 


dt 


d 


R{f,  f)  +  //■(/  X  g) 


—  a  ■ 


A 


+  p- 


X  g  -  /  X 


+ 


(m  X  /) 


d 


(m  X  /) 


(56) 


In  so  writing  (jsni)  we  have  deliberately  cast  it  into  a  form  that  eases  comparison  with  (El- 
In  the  appendix  we  set  up  an  apparatus  from  which  the  partial  derivatives  of  the  inertial 
terms  with  respect  to  the  orbital  elements  may  be  obtained.  We  also  show  that  some  of  these 
derivatives  vanish.  However,  a  complete  evaluation  of  the  inertial  input  to  the  planetary 
equations  in  the  ordinary  Lagrange  gauge  involves  a  long  and  tedious  calculation  which  we 
do  not  carry  out. 


3.5  Comparison  of  the  two  gauges 

One  of  the  powers  of  gauge  freedom  lies  in  the  availability  of  gauge  choices  that  simplify  the 
planetary  equations,  as  we  can  see  from  contrasting  (El  with  (El-  While  the  latter  equation 
is  written  under  the  customary  Lagrange  constraint  (i.e.,  for  elements  osculating  in  the  frame 
where  they  are  defined),  the  former  equation  is  written  under  a  nontrivial  constraint  called 
the  ’’generalised  Lagrange  gauge.”  The  simplicity  of  El  speaks  for  itself. 

By  identifying  the  parameters  Ci  with  the  Delaunay  variables,  one  arrives  from  El  and 
(El  to  the  appropriate  Delaunay-type  equations  (see  Appendix  I  to  Efroimsky  &  Goldreich 
2003a).  The  Lagrange  equations  corresponding  to  and  to  may  be  derived  from 
each  of  these  two  equations  by  choosing  C  as  the  Kepler  elements  and  using  the  appro¬ 
priate  Lagrange  brackets.  These  Lagrange  equations  are  written  down  in  the  Appendix  to 
Efroimsky  &  Goldreich  (2003b)  to  which  we  refer  the  interested  reader. 

Although  the  planetary  equations  are  much  simpler  in  the  generalised  Lagrange  gauge 
than  in  the  ordinary  Lagrange  gauge,  some  of  these  differences  are  less  important  than 
others.  In  many  physical  situations,  though  not  always,  the  p  ^  and  p  terms  in  (El  are  of 
a  higher  order  of  smallness  compared  to  those  linear  in  p ,  and  therefore  may  be  neglected, 
at  least  for  sufficiently  short  times.® 


4  Planetary  Equations  and  Gauges  in  the  Hamilton- 
Jacobi  Approach 

In  this  section  we  demonstrate  that  the  derivation  of  planetary  equations  in  the  N-particle 
{N  >  3)  case,  performed  through  the  medium  of  Hamilton- Jacobi  method,  implicitly  con¬ 
tains  a  gauge-hxing  condition  not  visible  to  the  naked  eye.  We  present  a  squeezed  account  of 
our  study;  a  comprehensive  description  containing  technical  details  may  be  found  in  Efroim¬ 
sky  &  Goldreich  (2003a). 

The  Hamilton- Jacobi  analysis  rests  on  the  availability  of  different  canonical  descrip¬ 
tions  of  the  same  physical  process.  Any  two  such  descriptions,  ( <?,  p,  H{q,p))  and 

®  As  an  example  of  an  exception  to  this  rule,  we  mention  Venus  whose  wobble  is  considerable.  This  means 
that,  for  example,  the  ft  term  cannot  be  neglected  in  computations  of  circumvenusian  orbits. 
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{Q,  P,  H*{Q,P))  ,  correspond  to  different  parametrisations  of  the  same  phase  flow,  and 
both  obey  Hamilton’s  equations.  Due  to  the  latter  circumstance  the  infinitesimally  small 
variations 

d9  =  pdq  —  Hdt  (57) 

and 

dO  =  PdQ  -  H*dt  (58) 

are  perfect  differentials,  and  so  is  their  difference 

-  dW  =  dO  -  dO  =  PdQ  -  pdq  -  {H*  -  H)  dt  .  (59) 


Here,  vectors  g,  p,  Q,  and  P  each  contain  N  components.  Given  a  phase  flow  parametrised 
by  a  set  (g,  p,  H{q,p,t) )  ,  it  is  always  useful  to  simplify  the  description  by  a  canonical 
transformation  to  a  new  set  {Q,  P,  H*{Q,  P,t) )  ,  with  the  new  Hamiltonian  H*  being 
constant  in  time.  Most  advantageous  are  transformations  that  nullify  the  new  Hamiltonian 
H*  ,  because  then  the  new  canonical  equations  render  the  variables  {Q,  P)  constant. 
A  powerful  method  of  generating  such  transformations  stems  from  (EHI)  being  a  perfect 
differential.  It  is  sufficient  to  consider  W  to  be  a  function  of  the  time  and  only  two  other 
canonical  variables,  for  example  g  and  Q  .  Then  (P|)  may  be  written  as 

dW  dW  dW 

-  -  ^dQ  -  ^dq  =  PdQ  -  pdq  +  {H  -  H*)  dt  (60) 

at  oQ  oq 


from  which  it  follows  that 

^  dW  dW 

~  ~  ~dQ  ’  ^  ~ 


H{q,  p,  t)  + 


dW 

~dt 


H%Q,  P,  t) 


(61) 


Inserting  the  second  equation  into  the  third  and  assuming  that  H*{Q,  P,  t)  is  simply  a 
constant,  we  get  the  famous  Jacobi  equation 


H 


^  dW  ^ 
.  dq  J 


+ 


dW 


=  H* 


(62) 


whose  solution  furnishes  the  transformation-generating  function  W.  The  elegant  power  of  the 
method  becomes  most  visible  if  the  constant  H*  is  set  to  zero.  Under  this  assumption  the  re¬ 
duced  two-body  problem  is  easily  resolved.  Starting  with  the  three  spherical  coordinates  and 
their  canonical  momenta  as  (g,  p)  ,  one  arrives  to  canonically  conjugate  constants  {Q,  P) 
that  coincide  with  the  Delaunay  elements  {Qi ,  Pi)  =  (T ,  —Mo)  ;  {Q2 ,  P2)  = 

(G,  -cn);  (Q3,  P3)  =  (H,  -  D)  . 

Extension  of  this  approach  to  the  N-particle  problem  begins  with  consideration  of  a  dis¬ 
turbed  2-body  setting.  The  number  of  degrees  of  freedom  is  still  the  same  (three  coordinates 
g  and  three  conjugate  momenta  p),  but  the  initial  Hamiltonian  is  perturbed: 


d{H  +  AH) 
dp 


p 


d{H  +  AH) 
dq 


(63) 


While  in  PI)  -  (El  one  begins  with  the  initial  Hamiltonian  H  and  ends  up  with  H*  =  0  , 
the  method  may  be  extended  to  the  perturbed  setting  by  accepting  that  now  we  start 
with  a  disturbed  initial  Hamiltonian  H  +  AH  and  arrive,  through  the  same  canonical 
transformation,  to  an  equally  disturbed  eventual  Hamiltonian  H*  +  AH  =  AH  .  Plugging 
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these  new  Hamiltonians  into  PI|)  leads  to  cancellation  of  the  disturbance  AH  on  the 
right-hand  side,  whereafter  one  arrives  to  the  same  equation  for  W{q,  Q,  t)  as  in  the 
unperturbed  case.  Now,  however,  the  new  canonical  variables  are  no  longer  conserved  but 
obey  the  dynamical  equations 


dAH  ■  _  dAH 


(64) 


Because  the  same  generating  function  is  used  in  the  perturbed  and  unperturbed  cases,  the 
new,  perturbed,  solution  {q,  p)  is  expressed  through  the  perturbed  ’’constants”  Q{t)  and 
P{t)  in  the  same  manner  as  the  old,  undisturbed,  q  and  p  were  expressed  through  the  old 
constants  Q  and  P  .  This  form-invariance  provides  the  key  to  the  N-particle  problem:  one 
should  choose  the  transformation-generating  function  W  to  be  additive  over  the  particles 
and  repeat  this  procedure  for  each  of  the  bodies,  separately. 

Armed  with  this  preparation,  we  can  proceed  to  uncover  the  implicit  gauge  choice  made 
in  using  the  Hamilton-Jacobi  method  to  derive  evolution  equations  for  the  orbital  elements. 
To  do  this  we  substitute  the  equalities 


and 


Q  = 


P  =  - 


dAH 

dAH  dq 

dAH  dp 

dP 

dq  dP  ' 

dp  dP 

dAH 

dAH  dq 

dAH  dp 

dQ 

dq  dQ 

dp  dP 

into  the  expression  for  the  velocity 


dq 


.  _  .  dq  ■  dq  ■ 

-  m  ^  +  ap^  ■ 


This  leads  to 


Q  = 


dq 

m 


+ 


dq  dq  dq  dq\  d  AH 


dQ  dP  dP  dQ 


fdAH\ 


9  = 


q,t 


dq 

dq 

m 


+ 


dq  dp 
dQ  ^ 


dq  dp\  d  AH 

dP  dQ  J  dp 


(65) 

(66) 

(67) 

(68) 


where  we  have  taken  into  account  that  the  Jacobian  of  the  canonical  transformation  is  unity: 


dq  dp  dq  dp 


(69) 


To  establish  the  link  between  the  regular  VOP  method  and  the  canonical  treatment,  compare 
(El  with  (TO-  We  see  that  the  symplectic  description  necessarily  imposes  a  particular  gauge 
$  =  5  AH / dp  . 

It  can  be  easily  demonstrated  that  this  special  gauge  coincides  with  the  generalised 
Lagrange  gauge  El  discussed  in  subsection  2.2.  To  that  end  one  has  to  compare  the 
Hamilton  equation  for  the  perturbed  Hamiltonian  dH, 


Q 


d  {H  +  AH) 
dp 


p  + 


dAH 

dp 


(70) 


16 


with  the  definition  of  momentum  from  the  Lagrangian  (CZD, 


p  = 


d  {L{q,  q,  t)  +  AL{q,  q,  t)  ) 
dq 


=  q  + 


OAL 

dq 


(71) 


Equating  the  above  two  expressions  immediately  yields 


$  = 


\  dp  J 


q,t 


.  9<1  / 


(72) 


q,t 


which  coincides  with  (EDI)-  Thus  the  transformation  generated  by  W{q,  Q,  t)  is  canonical 
only  if  the  physical  velocity  q  is  split  in  a  fashion  prescribed  by  (IHHll.  i.e.,  if  (17^  is  fulfilled. 
This  is  exactly  what  our  Theorem  from  subsection  2.2  says. 

To  summarise,  the  generalised  Lagrange  constraint,  $  =  —  dAL/dq  ,  is  tacitly  in¬ 
stilled  into  the  Hamilton-Jacobi  method.  Simply  by  employing  this  method  (at  least,  in 
its  straightforward  form),  we  automatically  fix  the  gauge. ^  By  sticking  to  the  Hamiltonian 
description  we  sacrifice  gauge  freedom. 

Above,  in  subsection  2.3,  we  established  that  in  the  generalised  Lagrange  gauge  the 
momentum  coincides  with  g  .  We  now  can  get  to  the  same  conclusion  from  (P|).  (ED)  and 

(1Z2D; 


_  9  {L{q,  q,  t)  +  AL{q,  q,  t)) 
^  ~  dq 


=  q  -  ^  =  g 


(73) 


Thus,  implementation  of  the  Hamilton-Jacobi  theory  in  celestial  mechanics  demands  the 
orbital  elements  to  osculate  in  phase  space.  Naturally,  this  demand  reduces  to  that  of 
regular  osculation  in  the  simple  case  of  velocity-independent  AL  . 


5  Conclusions 

In  the  article  thus  far  we  have  studied  the  topic  recently  raised  in  the  literature:  the  planetary 
equations’  internal  symmetry  that  stems  from  the  freedom  of  supplementary  condition’s 
choice.  The  necessity  of  making  such  a  choice  constrains  the  trajectory  to  a  9-dimensional 
submanifold  of  the  12-dimensional  space  spanned  by  the  orbital  elements  and  their  time 
derivatives.  Similarly  to  the  field  theory,  the  choice  of  the  constraint  (=  the  choice  of  gauge) 
is  vastly  ambiguous  and  reveals  a  hidden  symmetry  instilled  in  the  description  of  the  N-body 
problem  in  the  language  of  orbital  elements. 

We  addressed  the  issue  of  internal  freedom  in  a  sufficiently  general  setting  where  a  per¬ 
turbation  to  the  two-body  problem  depends  not  only  upon  positions  but  also  upon  velocities. 
Such  situations  emerge  when  relativistic  corrections  to  Newton’s  law  are  taken  into  account 
or  when  the  VOP  method  is  employed  in  rotating  systems  of  reference. 

Just  as  a  choice  of  an  appropriate  gauge  simplifies  solution  of  the  equations  of  motion 
in  electrodynamics,  an  alternative  (to  that  of  Lagrange)  choice  of  gauge  in  the  celestial  me¬ 
chanics  can  simplify  orbit  calculations.  We  provided  one  such  example, a  satellite  orbiting  a 
processing  planet.  In  this  example,  the  choice  of  the  generalised  Lagrange  gauge  considerably 
simplifies  matters. 

^An  explanation  of  this  phenomenon  from  a  different  viewpoint  is  offered  in  Section  6  of  Efroimsky  (2003) 
where  the  Delaunay  equations  are  derived  also  through  a  direct  change  of  variables.  It  turns  out  that  the 
outcome  retains  the  symplectic  form  only  if  an  extra  constraint  is  imposed  by  hand. 
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We  have  explained  where  the  Lagrange  constraint  tacitly  enters  the  Hamilton-Jacobi 
derivation  of  the  Delaunay  equations.  This  constraint  turns  out  to  be  an  inseparable  (though 
not  easily  visible)  part  of  the  method:  in  the  case  of  momentum-independent  disturbances, 
the  N-body  generalisation  of  the  2-body  Hamilton-Jacobi  technique  is  legitimate  only  if  we 
use  orbital  elements  that  are  osculating.  In  the  situation  where  the  disturbance  depends 
not  only  upon  positions  but  also  upon  velocities,  another  constraint  (which  we  call  the 
’’generalised  Lagrange  constraint”)  turns  out  to  be  stiffly  embedded  into  the  Hamilton-Jacobi 
development  of  the  problem. 

Unless  a  specihc  constraint  (gauge)  is  imposed  by  hand,  the  planetary  equations  assume 
their  general,  gauge-invariant,  form.  In  the  case  of  a  velocity-independent  disturbance,  any 
gauge  different  from  that  of  Lagrange  drives  the  Delaunay  system  away  from  its  symplectic 
form.  If  we  permit  the  disturbing  force  to  depend  also  upon  velocities,  the  Delaunay  equa¬ 
tions  retain  their  canonicity  only  in  the  generalised  Lagrange  gauge.  Interestingly,  in  this 
special  gauge  the  instantaneous  ellipses  (hyperbolae)  osculate  in  phase  space. 

Briefly  speaking,  N-body  celestial  mechanics,  expressed  in  terms  of  orbital  elements,  is 
a  gauge  theory  but  it  is  not  strictly  canonical.  It  becomes  canonical  in  the  generalised  La¬ 
grange  gauge. 
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Appendix 


In  this  appendix  we  set  up  an  apparatus  from  which  one  may  evaluate  the  partial  deriva¬ 
tives  with  respect  to  the  orbital  elements  of  inertial  terms  that  appear  in  the  planetary 
equations  derived  in  sections  3.3  and  3.4.  We  then  show  that  some  of  these  derivatives  van¬ 
ish.  Following  that,  we  derive  explicit  expressions  for  each  derivative  of  R-  (/  x  ,  which 
provides  a  complete  analytic  evaluation  of  the  rotational  input  in  the  generalised  Lagrange 
gauge.  The  topic  is  further  developed  (and  the  appropriate  generalised  Lagrange  system  of 
equations  is  presented)  in  (Efroimsky  &  Goldreich  2003b). 

To  hnd  the  explicit  form  of  the  dependence  /  =  /(Gj,  t)  ,  it  is  conventional  to  introduce 
an  auxiliary  set  of  Gartesian  coordinates  q  ,  with  an  origin  at  the  gravitating  centre,  and 
with  the  first  two  axes  located  in  the  plane  of  orbit.  The  q  coordinates  are  easy  to  express 
through  the  major  semiaxis  a ,  the  eccentricity  e  and  the  eccentric  anomaly  E  : 

qi  =  a{cosE  —  e)  ,  q2  =  aVi  —  G  sinE  ,  qs  =  0  ,  (74) 

where  E  itself  is  a  function  of  the  major  semiaxis  a  ,  the  eccentricity  e  ,  the  mean  anomaly 
at  epoch,  Mq  ,  and  the  time,  t .  The  time  dependence  is  realised  through  the  Kepler  equation 

E  —  e  sinE  =  M  ,  (75) 

where 

M  =  Mo  +  f'  dt  .  (76) 
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The  inertial-frame-related  position  of  the  body  reads: 

r  =  /  (fl,  i,  u,  a,  e,  Mo  ;  t)  =  R(r2,  z,  a;)  q(a,  e,  E{a,  e,  Mq,  f) )  ,  (77) 

R(n,  z,  a;)  being  the  matrix  of  rotation  from  the  orbital-plane-related  coordinate  system 
q  to  the  fiducial  frame  {x,  y,  z)  in  which  the  vector  r  is  defined.  This  rotation  is 
parametrised  by  the  three  Euler  angles:  inclination,  i  ;  the  longitude  of  the  node,  ;  and 
the  argument  of  the  pericentre,  to  . 

In  the  unperturbed  two-body  setting  the  velocity  is  expressed  by 


d  ( dE  \  ■'  ( do '' 

g  =  f  i,  uj,  a,  e,  Mo  ]  t)  =  f  —  J  R(f],  i,  u)  i  — 

V  /  a,  e,  Mo  \  Z 

One  can  similarly  calculate  partial  derivatives  of  /  with  respect  to  Mo  : 


d  (  dE 

f  (O,  i,  uj,  a,  e.  Mo]  t)  = 


dM, 


dMo 


R(SJ,  !,!..)  (|| 


^  /  a,e,t 

whence  it  becomes  evident  that  df/dMo  is  parallel  to  g  and,  hence, 

df 


g  X 


dMo 


=  0  . 


i,  LO,  a,  e,  t 


(78) 


(79) 


(80) 


By  a  similar  trick  it  is  possible  to  demonstrate  that  d{f  x  g)/dMo  is  proportional  to 
d{f  X  g)/dE  and,  therefore,  to  d{f  x  g)/dt  .  Hence,  this  derivative  vanishes  (because  in 
the  two-particle  case  the  cross  product  /  x  g  is  an  integral  of  motion).  This  vanishing  of 
d{f  X  g)/dMo  ,  along  with  (jHOI),  implies: 


0 


(81) 


In  the  situation  when  the  parameters  Oj  are  implemented  by  the  Delaunay  elements,  a 
similar  sequence  of  calculations  leads  to 


g  X 


n,  Lj,  L,  G,  t 


0 


(82) 


and,  appropriately,  to 


\dMo)  n,uj,L,G,H,t 


0  . 


(83) 


We  can  proceed  much  farther  in  the  generalised  Lagrange  gauge,  at  least  in  so  far  as 
derivatives  of  the  rotational  input  J//i  =  f  x  g  are  concerned.  (We  remind  that  here  and 
everywhere  y  stands  for  the  reduced  mass,  while  p,  denotes  the  precession  rate.) 

As  we  proved  above,  this  cross  product  is  independent  of  Mo  and,  hence. 


^  djf  X  g) 
^  ■  dMo 


=  0  . 


(84) 


19 


Since  J  is  orthogonal  to  the  orbit  plane,  it  is  invariant  under  rotations  of  the  orbit 
within  its  plane,  whence 

^  d{f  X  g)  .  . 


To  continue,  we  note  that  in  the  two-body  setting  the  ratio  J//i  ,  is  known  to  be  equal  to 
w  where  w  is  a  unit  vector  perpendicular  to  the  unperturbed  orbit’s  plane. 
Moreover,  in  planet-associated  noninertial  axes  {x,  y,  z)  with  corresponding  unit  vectors 
(x,  y,  z) ,  the  normal  to  the  orbit  is  expressed  by 

w  =  X  sin  i  sin  hi  —  y  sin  i  cos  hi  -|-  z  cos  i  .  (86) 


Hence, 


d  1 

#  t  .  _ 

[fyg] 

1  ^  ^  d  \ 

_  It  .  WT  _ 

(^a(l  -  e2)  ^ 

ll  -  e2 

da 

-  IjL  W 

da 

2  V 

a 

^  9  1 

§  1  .  _ 

[f  x  a) 

^  d  \ 

-  I  1  .  WT  _ 

[^a{l  -  e2)  ) 

1  a 

de 

-  Li/  W 

de 

■  =  -v 

1  -  e2 

where  sin  z  sin  hi  —  fiy  sin  i  cos  fl  +  yz  cos  i  is  the  orthogonal-to-orbit 

component  of  the  precession  rate.  The  remaining  two  derivatives  look: 


a  (/ 


' a  {1  —  e^)  {yx  sin  z  cosh!  -|-  yy  sin  z  sin  hi} 


and 


-  a  (/ X  9) 


a  (1  —  e^)  {yx  cos  z  sin  hi  —  yy  cos  z  cosh!  —  yz  sin  z} 


(90) 


As  for  the  derivatives  of  a  ■  /  ,  they  may  be  calculated  directly  from  the  expression  for 
ctz,  z,  a,  e.  Mo]  t)  presented  above.  However,  the  resulting  expressions  are  cumbersome 
so  we  do  not  present  them  here. 
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